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Oscillating magnetocaloric effect of a 2D non-relativistic diamagnetic material 
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Among the magnetic materials, those with ferromagnetic character are, by far, the most studied 
in what concerns applications of the magnetocaloric effect. However, recently, diamagnetic materi¬ 
als received due attention never received before, and an oscillatory behavior, analogous to the de 
Haas-van Alphen effect, has been found. The present effort describes in details the magnetocalo¬ 
ric properties of a 2D non-relativistic material (a Gold thin film, for instance), where oscillations, 
depending on the reciprocal magnetic field 1 /B, are found. A comparison of the magnetic entropy 
change per electron for some cases is presented and we found « HR 1 kn (@109.3 K) for graphenes, 

« 10~ 5 kn (@0.7 K) for 2D Gold and « 1CF' ku (@0.7 K) for 3D Gold. 
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I. INTRODUCTION 

The magnetocaloric effect is an interesting property in 
which magnetic materials, under a magnetic field change, 
are able to exchange heat with a thermal reservoir (con¬ 
sidering an isothermal process), or even change its tem¬ 
perature (considering an adiabatic process). This ef¬ 
fect is completely analogous to the compression-expasion 
thermal-mechanical cycle; and therefore the main pur¬ 
pose of the scientific community is to built a thermo- 
magnetic machine, to substitute, in a near future, those 
standard, non-economical and non-environmental frien¬ 
dly Freon-like refrigerators [l|. 

Thus, the scientific community dealing with caloric ef¬ 
fects (baro-, magneto- and electro-), focused attention 
to applied purposes and therefore, to this end, ferro- 
materials are the most useful, since near phase transi¬ 
tions the caloric effects are maximized[l|. However, re¬ 
cently, the magnetic entropy change AS(T, AB) and the 
adiabatic temperature change AT(T, AB) of 3D non- 
relativistic diamagnetic materials have been described in 
details^, 0|- An oscillatory behavior was found, due to 
the crossing of the Landau levels with the Fermi level, 
analogously to the de Haas-van Alphen effect Q. This 
new (oscillating) magnetocaloric effect has interesting fe¬ 
atures and, as suggested, can be used as sensible magne¬ 
tic field sensor Q- 

However, in spite to be interesting, this oscillating ef¬ 
fect occurs at high values of magnetic field (few or more 
Teslas) and low temperatures (c.a. 1 K). In order to 
further understand and optimize the effect, a relativistic 
2D material (a graphene) was studiedjsj; and, surprisin¬ 
gly, a very similar oscillations were found, however, at 
very comfortable values of temperature: c.a. 100 K (and 
few Teslas). To deeper understand the magnetocaloric 
properties of graphenes, more recently the influence of 
a longitudinal electric field was described Q and verified 


that the caloric potentials can be easily rules by this elec¬ 
tric field. The electrocaloric effect of graphens was also 

explored Q- 

What does happen for a material that has low dimen¬ 
sion (2D), like graphene, but is non-relativisit in nature, 
as the first case above described? This is the aim of the 
present effort, to present the magnetic entropy change of 
a 2D non-relativistic film. To this purpose, we derive the 
density of states (details are in the appendix), and then, 
from known thermodynamic relationships, the magnetic 
entropy change is obtained. 


II. MAGNETIC ENTROPY CHANGE 

The magnetic entropy S(T,B), as a function of tem¬ 
perature T and magnetic field B , depends on the grand 
canonical partition function Z(T, B) through: 

S(T,B) = k B -^[TlnZ(T,B)] (1) 


where 


POO 

In Z(T,B)= / de gi(e) In (l + ze~^ e ^ , 

J o 


( 2 ) 


2 = e^ 3 is the fugacity, g is the chemical potential and 
/3 = l/k B T. In addition, gi(e) is the one-particle density 
of states, given by (see Appendix): 

5i(e) = So + Ss(e) (3) 


where 


L 2 2m 
90 = A^H 2 


(4) 


is the zero-field density of state (with no spin degene¬ 
racy), and 


+oo 

g B (s) = 2g 0 ^2(-l) 1 cos 
1=1 



( 5 ) 
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From equation [3J it is easy to see that the logarithm 
of the grand partition function has two contributions: 

lnZ(T,B) =lnZ 0 (T)+lnZ B (T) ( 6 ) 

The first one gives the zero-field entropy and, as will be 
shown further in this text, does not contribute to the 
magnetic entropy change. Thus, let us focus our atten¬ 
tion to the second term that, from equation [5j is: 


roo +°° / j \ 

lnZs(T) = J de2g 0 ^(-1)* cos ^ e) In (l + ze~ 0e ) 


(7) 


l=i 

After integration by parts (two times), we achieve: 

In Zb (T) = 2 5 o£(- 1 )* {g( £ )| 0 °° + S(e)| 0 °° 


1 UbBP V 

A\ In ) 


3{T,B)\ (8) 


where 


G(£) = ^ sin (_^_ e ) ln(l + ze-*) 


(9) 


S(£) = 


HbB f In \ 1 

In \p B B J z 1 e@ e + l 


( 10 ) 


and 


3(T,B) = 


r cos (^g £ ) . 

Jo cosh 2 (fa)) 

e F » fcB T 47 r 2 l C0S (i7TB^) 

~^ VbB/ 3 2 sinh {_w\ 

\ii b B0 J 


( 11 ) 


Above, the condition £p k B T implies to p, ~ ep and 
therefore z 1. As a consequence, G(oo), G(0) and 
9 (oo) are zero; and 


s(0) - 


VbB 

In 


( 12 ) 


Equation |7] can then be written with two contributions: 

In Z b (T) = InZ n B °(T) + In Z° B (T) (13) 

where 


\nZ n B °{T) = -\g 0 ^ B B) 2 
6 


(14) 


has a non-oscillatory character and 

+°° c_iy 

In Z° B (T) = —2gag, B B — 7 — cos 


1=1 


l 


f- 


In 


\PbB 


£ F 


1 


sinh(a:;) 

(15) 


has an oscillating behavior. Above, Xi = lx and 


x = n 


2 


k B T 

g B B 


Thus, equation [ 6 ] reads as: 


(16) 


lnZ(T,B) =\nZ 0 (T) + hiZ% o (T)+lnZ o B (T) (17) 

where the second and third terms are given by equations 
M and [15j respectively. From equation [TJ it is possible 
to obtain the entropy of the system, that reads as: 

S(T,B) = S 0 (T) + S° b (T) (18) 


since S B ° = 0. However, we are interested to obtain the 
magnetic entropy change A S(T, A B) = S(T, B) — S(T , 0) 
that, on its turn, is: 


A S(T, A B) 

= aSd(T) = — 2NkB— costnn)*! (x) 
n 

(19) 

since S° B=Q {T) 

= 0. Above, 



£F 

n =- 

( 20 ) 


PbB 

and, for convenience, we re-write x = n 2 tn , where 



t- kBT 

( 21 ) 


£F 

In addition, 

7( X ) - xL( - x) 
sinh(a;) 

( 22 ) 

and L(x) is the Langevin function: 



L(x) = coth(x) — — 

X 

(23) 


As usually considered for thermodynamic quantities with 
this kind of oscillations [§, jjj, we considered only 1 = 1 
term: the hyperbolic sine function at the denominator of 
equation [22] rapidly damping the summation. Of course, 
it is an approximation; however, after this assumption, 
we can understand the physics behind this model. 

Let us consider a Gold thin film, with Fermi energy 
Ef = 3.62 eV. From equation [5D] it is possible to see 
that B{n = 1) = 6.2522 x 10 4 T and therefore com¬ 
pletely out of a laboratory range. On the other hand, 
B{n = 10 4 ) = 6.2522 T, and therefore this is the or¬ 
der of magnitude of n we must consider. Note the 
magnetic field change needed to invert the magnetic en¬ 
tropy change (from normal/negative to inverse/positive) 
is B[n = 10 4 + 1) = 6.2515 T, i.e., a difference of 0.7 mT. 
The temperature dependence of the magnetic entropy 
change lies on the 7{x) function and this behavior is pre¬ 
sented in figure [I] This function peaks in x = n 2 tn =1.6 
and therefore, considering n = 10 4 , the maximum mag¬ 
netic entropy change occurs at T = 0.7 K. 

The dependence of the magnetic entropy change as a 
function of n is presented in figure^ for that temperature 
that maximizes this effect (considering n = 10 4 ), i.e., 0.7 
K. It has an oscillatory behavior due to the cosine term 
on equation [Hi] (note the envelope of these oscillations is 
equation 1191 without the cosine term). 























3 


tion as these presented in this contribution, the magne- 



Figura 1. Oscillating magnetic entropy change per electron as 
a function of temperature, for a 2D Gold film. Note only 0.7 
mT of change on the magnetic field change is enough to invert 
the caloric effect from normal/negative to inverse/positive. 



ture dependence of the magnetic entropy change, both 
are driven by T (y). Since the argument of this function 
is (almost) the same for both cases, then T(x) peaks ap¬ 
proximately at the same value of temperature (0.7 K), for 
the 2D and 3D models. Thus, figure [I] is a good approxi¬ 
mation for the behavior of the magnetic entropy change 
as a function of temperature for the 3D non-relativistic 
case; except for the fact the 2D case is almost 2 orders of 
magnitude bigger than the 3D case (for n = 10 4 ) - this 
fact will be clearer further in the text. 

On the other hand, the magnetic entropy change as 
a function of n has important differences. In spite of 
the period of oscillations be (almost) the same for both 
cases (the differences lie on the Fermi energy values and 
the 1/4 factor-see equation [25] and [13). the envelope of 
these oscillations are remarkable different - the envelope 
is the entropy change without the cosine terms. Note the 
maximum amplitude for the 3D case rapidly goes to zero, 
while for the 2D case it keeps almost constant and then, 
for large values of n, suddenly goes to zero. However, 
as discussed above, n = 10 4 corresponds to B = 6.2522 
T, a range of magnetic field close to standard laboratory 
values; and therefore, for n = 10 4 , the magnetic entropy 
change, per electron, for a 2D Gold film is two orders of 
magnitude bigger than the 3D Gold case. 


Figura 2. Oscillating magnetic entropy change per electron as 
a function of n, inversely proportional to the magnetic field 
B. The oscillatory behavior is evident and has a remarkable 
difference comparing 2D and 3D models. The envelope rules 
the amplitude of the oscillations. For the sake of clearness of 
the figure, above n = 10 only the envelope is shown; and the 
oscillations are suppressed. 


III. COMPARISON WITH A 3D 
NON-RELATIVISTIC MATERIAL 

The magnetocaloric effect of diamagnetic materials be¬ 
gun to be discussed very recently, considering a 3D non- 
relativistic Gold specimen^. After a similar evalua- 


IV. COMPARISON WITH A 2D RELATIVISTIC 
MATERIAL: GRAPHENES 

More recently the caloric properties of graphenes have 
been described [5j-(8|; in a similar fashion as presented 
above. The magnetic entropy change per graphene area 
was evaluated and reads as[5j: 

AS(T, AB) = 2k B — cos(7rm)T(z) (27) 
m 

where No = 10 16 m -2 is the density of charge carriers^, 

0 , 


m = Nq 


B ’ 


(28) 
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Figura 3. Oscillating magnetic entropy change per electron 
for a graphene, as a function of temperature. 


(f >o = 7 xh/e = 2.06 x 10 15 Tm 2 is the magnetic flux 
quantum, 


m ksT 
N 0 v F 


(29) 


v F = hv F /2ny/Ngn = 9.43 x 10 _38 Jm 2 and v F = 10 6 
m/s is the Fermi velocity. Note both, to and x are di¬ 
mensionless. Interestingly, T (z) is also given by equation 

M 

The entropy change of graphenes (equation E71) . pro¬ 
vides a huge difference to what was discussed until this 
point. Note B(m = 1) = 20.6 T and B(m = 3) = 6.9 T, 
i.e., within the laboratory range. Thus, for this graphene 
comparison, we will focus to m ~ 1, instead of those 
n = 10 4 for a non-relativistic material. In addition, 
the temperature dependence of this entropy change for 
graphenes is also on T (z) function of equation l22l that 
peaks at z = 1.6 and therefore, for to = 1, the maximum 
magnetic entropy change occurs at T=109.3 K. Figure 
[3] summarizes these considerations and, comparing with 
figure [I] it is obvious the huge difference, mainly in the 
temperature in which the oscillating behavior occurs. 

Analogously to figure [21 the entropy change as a func¬ 
tion of to for graphenes is presented in figure @]-bottom, 
for values of magnetic field within the laboratory range. 
Note the oscillations rapidly goes to zero and few peri¬ 
ods of oscillations can be observed. The most interesting 
point is the temperature and magnetic held in which the 
effect is; i.e., completely reasonable in what concerns a 
standard laboratory. For the sake of comparison, the en¬ 
tropy change as a function of n (see equation 1191) is pre¬ 
sented in figure 01-top, for a 2D non-relativisitc material. 


These huge differences came from the relativistic 
energy spectra for the graphene, in which the Landau 
levels goes proportionally to y/B j, where j is the Lan¬ 
dau level index; while, as known and shown in Appendix, 




m = N A /B (dimensionless) 


Figura 4. Oscillating magnetic entropy change per electron for 
a (top) 2D non-relativistic film-Gold and (bottom) graphene, 
as a function of the inverse magnetic field. 


for a non-relativistic case, the energy spectra mimics an 
harmonic oscillator. 

As a final remark, at standard laboratory values of 
magnetic held change (c.a. 6 T), the entropy change per 
electron for a 2D relativistic material (graphene), is of 
the order of 10 _1 ke (@109.3 K), while for a 2D and 
3D non-relativistic material (Gold, for instance), it is, 
respectively, of the order of 10 -5 ke and 10“' ke (@0.7 
K). 


V. CONCLUSIONS 

In what concerns caloric effects, the scientihc com¬ 
munity has been focusing attention to ferro- materials, 
since these effects are maximized close to phase transiti¬ 
ons. However, recently, the magnetocaloric effect of a 3D 
non-relativistic diamagnetic materials has been descri¬ 
bed and, surprisingly, interesting and new features arose, 
as, for instance, an oscillatory behavior as a function of 
the reciprocal magnetic held 1 /B, analogously to the de 
Haas-van Alphen effect. This effect was further deve¬ 
loped considering a relativistic material: a 2D carbon 
sheet, i.e., a graphene; and a similar behavior was found, 
however, for more comfortable values of temperature (c.a. 
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109 K, in contrast to the low values -c.a. 1 K- found for 
the non-relativistic case). 

In a sequence of the results mentioned above, the pre¬ 
sent effort dealt with a non-relativistic 2D material, na¬ 
mely a Gold thin film. We found the characteristic os¬ 
cillatory behavior already mentioned for those other di¬ 
amagnetic materials (compare equations [TDJ HH and [271 , 
however, with important aspects that must be empha¬ 
sized: (i) the temperature dependence for those three 
different cases lies on the T(a;) function (see equation 
1221) - independently on the dimension and nature, i.e., 
either relativistic or non-relativistic, (ii) all cases have 
dependence on a cosine function, with a period proporti¬ 
onal to the reciprocal magnetic field 1/B and (iii) for 
those 2D cases the entropy change is proportional to 
B cos(l/-B)T(l/f?), while for the 3D case it depends on 
B 3 / 2 cos(l/-B)T(l/i?). From this, it is possible to ins¬ 
pect that the entropy change for a d-dimensional system 
is B a ! 2 cos(l/f3)T(l/B). This proposal needs further 
analysis and will be published elsewhere. 


Apendice A: Density of States 


The one-particle density of states 51(e) follows easier 
from the Laplace transformation of the canonical parti¬ 
tion function Z^(/3) in the Boltzmann limit: 

1 r/3+ioo 

9i(e) = — e^Z?(P')d/3' (Al) 

Zm Jp- ioo 

where, for non-interacting systems, a simple grand 
canonical-canonical partition functions relationship is: 

In Z(T,B) = zZ?(/3). (A2) 

Above, z = is the fugacity and /.4 the chemical po¬ 
tential. Thus, at this point, we need to obtain the grand 
partition function in the Boltzmann limit: 


In Z(T,B) = Y J ze~ 0£n (A3) 

n 


but before, the energy spectra e n for this model. 

The Hamiltonian of this model, i.e., a 2D non- 
relativistic electron gas with a transversal magnetic field 
B = Bk 7 can be written as: 


‘ K= 2 m^ + ^ Py + eBX ^ 


(A4) 


and then the Schrodinger equation reads as: 


h 2 d 2 1 

-^,d^ + 2 mUj2{x + Xo)2 ^x) = £ n<t>n(x) (A5) 


where 

hk v eB ( 1 

Xq = ——, to = — and e n = nw \ n + — 

eB m \ 2 

( A6 ) 

Above, k y = 2ijn y /L (where n y = 0,1, 2, • • •), is related 
to the translational symmetry along the y axis; and the 
energy spectra represents the Landau levels, where n is 
the Landau level index. Note the harmonic oscillator of 
equation IA5I has several centers xo, that depends on n y 
and B. Since these centers must be within the considered 
plan, i.e., 0 < Xq < L then 0 < n y < g, where 5 = 
I?eB/h is the degeneracy of the Landau level n (note 
for each n there are n y possibilities). 

Thus, equation [A3] can be rewritten as: 


OO 

In Z(T , B) = zg exp 

n=0 

L 2 eB 1 

— g - 

2 h sinh(;y) 


( n + 


(A7) 


where y = hbB(3. Thus, a simple comparison of equation 
IA7I and IX2l leads to: 




L 2 eB 1 
2 h sinh(j/) 


(AS) 


From the above and equation [AT] it is possible to write 
the density of states we are looking for: 


5i(e) 


L 2 eB 
2 h 


1 

27T* 




sinh(/i ibB(3') 



(A9) 


that reads as: 


5i(e) =50 + 9b (e) 


where 


L 2 2m 

90 = ^K 2 


(A10) 


(All) 


is the density of state (with no spin degeneracy), of the 
non-perturbed 2D electron gas; and 


+ OO 

9B (e) = 2 5 o X](~ 1) Z cos 
1 = 1 



(A12) 
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